Non-Abelian Chern-Simons models with discrete gauge groups on a lattice 
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We construct the local Hamiltonian description of the Chern-Simons theory with discrete non- 
Abelian gauge group on a lattice. We show that the theory is fully determined by the phase factors 
associated with gauge transformations and classify all possible non-equivalent phase factors. We also 
construct the gauge invariant electric field operators that move fiuxons around and create/anihilate 
them. We compute the resulting braiding properties of the fiuxons. We apply our general results to 
the simplest class of non-Abelian groups, dihedral groups Dn- 



I. INTRODUCTION. 

Realization of a quantum computer would be very im- 
portant for the solution of hard problems in classical 
mathematicsi and even more important for the solu- 
tion of many quantum problems in field theory, quantum 
chemistry, material physics, etc.^-'^''* However, despite a 
large effort of many groups the realization remains far 
away because of the conflicting requirements imposed by 
scalability and decoupling from the environment. This 
dichotomy can be naturally resolved by the error free 
computation in the topologically protected space^-^ if the 
physical systems realizing such spaces are identified and 
implemented. Though difficult, this seems more realistic 
than attempts to decouple simple physical systems imple- 
menting individual bits from the environment. Thus, the 
challenge of the error free quantum computation resulted 
in a surge of interest to physical systems and mathemat- 
ical models that were considered very exotic before. 

The main idea of the whole approach, due to Kitae\iS, 
is that the conditions of long decoherence rate and scal- 
ability can be in principle satisfied if elementary bits are 
represented by anyons, the particles that indergo non- 
trivial transformations when moved adiabatically around 
each other (braided)2ii2iii and that can not be distin- 
guished by local operators. One of the most famous ex- 
amples of such excitations is provided by the fractional 
Quantum Hall EffectiLi^. The difficult part is, of course, 
to identify a realistic physical system that has such exci- 
tations and allows their manipulations. To make it more 
manageable, this problem should be separated into dif- 
ferent layers. The bottom layer is the physical system 
itself, the second is the theoretical model that identifies 
the low energy processes, the third is the mathematical 
model that starts with the most relevant low energy de- 
grees of freedom and gives the properties of anyons while 
the fourth deals with construction of the set of rules on 
how to move the anyons in order to achieve a set of uni- 
versal quantum gates (further lies the layer of algorithms 
and so on). Ideally, the study of each layer should pro- 
vide the one below it a few simplest alternatives and the 
one above the resulting properties of the remaining low 
energy degrees of freedom. 



In this paper we focus on the third layer: we discuss 
a particular set of mathematical models that provides 
anyon excitations, namely the Chern Simons gauge the- 
ories with the discrete gauge groups. Generally, an ap- 
pealing realization of the anyons is provided by the fiuxes 
in non-Abelian gauge theories iA2*ii. The idea is to en- 
code individual bits in the value of fiuxes belonging to 
the same conjugacy class of the gauge group: such fiuxes 
can not be distinguished locally because they are trans- 
formed one into another by a global gauge transformation 
and would be even completely indistiguishable in the ab- 
sence of other fluxes in the system. Alternatively, one can 
protect anyons from the adverse effect of the local opera- 
tors by spreading the fiuxes over a large region of space. 
In this case one does not need to employ a non-Abelian 
group: the individual bits can be even encoded by the 
presence/absence of fiux in a given area, for instance 
a large hole in the lattice. Such model3iSiiiiiiiS4i& are 
much easier to implement in solid state devices but they 
do not provide a large set of manipulation that would 
be sufficient for universal quantum computation. Thus, 
these models provide a perfect quantum memory but not 
a quantum processor. On the other hand, the difficulty 
with the fiux representation is that universal computa- 
tion can be achieved only by large non-Abelian groups 
(such as Ac,) that are rather difficult to implement, or if 
one adds charge motion to the allowed set of manipula- 
tions. Because the charges form a non-trivial represen- 
tation of the local gauge group, it is difficult to protect 
their coherence in the physical system which makes the 
latter alternative also difficult to realize in physical sys- 
tems. The last alternative is to realize a Chern-Simons 
model where on the top of the conjugacy transformations 
characteristic of non-Abelian theories, the fiuxes acquire 
non-trivial phase factors when moved around each other. 
We explore this possibility in this paper. 

Chern-Simons theories with finite non-Abelian gauge 
groups have been extensively studied for a continuous 
2-1-1 dimensional space-time. Unlike continous gauge- 
group, discrete groups on a continous space allow non- 
trivial fluxes only around cycles which cannot be con- 
tracted to a single point. So in a path-integral approach, 
and by contrast to the case of continuous groups, the 
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integration domain is restricted to gauge-field configura- 
tions for which the local flux density vanishes everywhere. 
Such path integrals were introduced, analyzed, and clas- 
sified in the original paper by Dijkgraaf and Witter^''. 
They showed that for a given finite gauge group G, all 
possible Chern-Simons actions in 2 + 1 dimensions are 
in one to one correspondence with elements in the third 
cohomology group H^{G, U{1)) of G with coefficients in 
C/(l). They also provided a description in terms of a 
2-1-1 lattice gauge theory, where space-time is tiled by 
tetrahedra whose internal states are described by just 
three independent group elements {g, h, k) because fluxes 
through all triangular plaquettes vanish. Elements of 
H^{G,U{1)) are then identified with functions a{h,k,l) 
that play the role of the elementary action for a tetra- 
hedron. This description turns out to be rather cumber- 
some because q;(/i, fc, I) does not have specific symmetry 
properties, so that the definition of the total action re- 
quires to choose an ordering for all the lattice sites, which 
cannot be done in a natural way. As a result, it seems 
very difficult to take the limit of a continuous time that 
is necessary for our purposes because physical implemen- 
tations always require an explicit Hamiltonian form. 

In principle, the knowledge of an elementary action a 
in H^{G, U{1)) allows to derive all braiding properties of 
anyonic excitations (i.e. local charges and fluxes). This 
has been done directly from the original path-integral 
formulationSLSS, using as an intermediate step the rep- 
resentation theory of the so-called quasi-quantum double 
associated to the group G^^. This mathematical struture 
is completely defined by the group G and a choice of any 
element a in H^{G,U{1)). Using this later formalism, 
which bypasses the need to construct microscopic Hamil- 
tonians, many detailed descriptions of anyonic properties 
for various finite groups, such as direct products of cyclic 
groups, or dihedral groups, have been given by de Wild 
Propitius^'*. Unfortunately, these general results do not 
provide a recipe how to contruct a microscopic Hamilto- 
nian with prescribed braiding properties. 

To summarize, in spite of a rather large body of the- 
oretical knowledge, we still do not know which Chern- 
Simons models with a finite gauge group can in principle 
be realized in a physical system with a local Hamiltonian, 
that is which one can appear as low energy theory. To 
answer this question is the purpose of the present paper. 

The main ingredients of our construction are the fol- 
lowing. We shall be mostly interested in a Hilbert space 
spanned by dilute classical configurations of fluxes be- 
cause it is only these configurations that are relevant for 
quantum computation that involves flux braiding and fu- 
sion. Furthermore, we expect that a finite spacial sep- 
aration between fluxons may facilitate various manipu- 
lations and measurements. Notice in this respect that 
in theories with discrete group symmetry in a contin- 
uous space^SiSiiS^j non-trivial fluxes appear only thanks 
to a non-trivial topology of the ambiant space^^ and thus 
are restricted to large well separated holes in a flat 2D 
space. The second feature of our construction is that 



gauge generators are modified by phase-factors which 
depend on the actual flux configuration in a local way. 
This is a natural generalization of the procedure we have 
used recently for the construction of Chern-Simons mod- 
els with Z N symmetry groupi^iSI. Note that in the most 
general situations, the representation of the local gauge 
group in a discrete Chern-Simons theory becomes projec- 
^j^g ^0,2ij22j23j24 ^ This would bc inappropriate for a ro- 
bust implementation because projective representations 
lead to degenerate multiplets that are strongly coupled to 
local operators, and therefore become very fragile in the 
presence of external noise. We shall therefore restrict our- 
selves to the models were no projective representations 
occur. In practice, they are associated with the non- 
trivial elements of the group i7^(G, [/(!)), and it turns 
out that for some interesting classes of groups such as the 
dihedral groups Dn with N odd, i?^(G, [/(!)) is triviaPi, 
so this restriction is not too important. As we show in 
Section lTnl these assumptions allow us to find all the pos- 
sible phase factors associated with the gauge transforma- 
tions in terms of homomorphisms from the subgroups of 
G that leave invariant a fixed element of G under conju- 
gacy into U{1). The last step is to construct a set of local 
gauge-invariant operators corresponding to the following 
elementary processes: moving a single fluxon, creating 
or annihilating a pair of fluxons with vacuum quantum 
numbers, and branching processes were two nearby flux- 
ons merge into a single one (or conversely). We shall see 
that the requirement of local gauge invariance leads to a 
relatively mild constraint for the possible phase-factors, 
leaving a rather large set of solutions. 

The main result of this work is twofold. First, we pro- 
vide an explicit computation of holonomy properties of 
fluxes in a set a models based on dihedral groups. Of spe- 
cial interest is the simplest of them, D3, which is also the 
permutation group of 3 elements S3 . This group belongs 
to a class which is in principle capable of universal quan- 
tum computatiooii. This part is therefore connected to 
the upper layer (designing a set of univeral gates from 
properties of anyons) in the classiflcation outlined above. 
But our construction of a local Hamiltonian version for 
a Chern-Simons model on a lattice provides some guide- 
lines for possible desirable physical implementations. 

The plan of the paper is the following. Section II 
is mostly pedagogical, providing the motivation for our 
general construction through the simplest example of 
a Chern-Simons gauge theory, namely the non-compact 
U{1) model. In Section III we formulate general condi- 
tions on local Chern-Simons phase factors that satisfy 
gauge invariance conditions. In Section IV we discuss 
the construction of the electric field operator and derive 
condition on the phase factor that allows one a gauge 
invariant fluxon dynamics. In Section V we discuss the 
fluxon braiding properties and derive the Chern-Simons 
phase factors associated with the braiding. In Section VI 
we apply our results to the simplest non-Abelian groups 
Dn. Finally, Section VII gives conclusions. Some techni- 
cal points relevant for the general construction have been 
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relegated to Appendix A, and Appendix B discusses some 
interesting features of the torus geometry. Although this 
geometry is not easy to implement in experiments, it is 
conceptually very interesting since it is the simplest ex- 
ample of a two-dimensional space without boundary and 
with topologically non-trivial closed loops. 

II. OVERVIEW ON ABELIAN CHERN-SIMONS 
MODELS 

To motivate the construction of the present paper, it 
is useful to discuss some important features of Chern- 
Simons gauge theories in the simpler case of Abelian 
groups. For this purpose, we shall consider here as an 
illustration the model based on the continuous Abelian 
group with one generator in its non-compact version. Of 
course, our main purpose here is to address finite groups, 
but as we shall discuss, this non-compact U{1) model 
contains already the key ingredients. On a 2 -f 1 dimen- 
sional space-time, it is defined from the following La- 
grangian density: 

C = ^X{Al + AD - ^^iB' + v{A.,Ay - AyA,,) (1) 

where B = dxAy — dyA^ is the local magnetic field (a 
pseudo-scalar in 2 -I- 1 dimensions) and dots stand for 
time-derivatives. We have used the gauge in which the 
time component Aq of the vector potential is zero. Be- 
cause of this, we shall consider only invariance under 
time-independent gauge transformations in this discus- 
sion. These are defined as usual by Ap ^ Ap + dpf , 
where f{x,y) is any time-independent smooth scalar 
function of spacial position. Under such a gauge trans- 
formation, the action associated to the system evolution 
on a two-dimensional space manifold M during the time 
interval [ti,i2] varies by the amount AS where: 

Because / is time-independent, the integrand 
is a total time-derivative, so we may write 
AS = iy{I{A2j) - liAiJ)), where A, denotes the 
field configuration at time ti, i ^ 1,2, and: 

In the case where M has no boundary (and in particular 
no hole), we may integrate by parts and get: 

When ^ 0, this modifies in a non-trivial way the 
behavior of the corresponding quantized model under 
time- independent gauge transformations. One way to 
see this is to consider the time-evolution of the system's 



wave- functional ^{A, t). In a path-integral approach, the 
probability amplitude to go from an initial configuration 
Ai{r) at time ti to a final A2{r) at time t2 is given by: 

A21 = J VA{r,t) exp{-^5(A)} 

where fields A{r, t) are required to satisfy boundary con- 
ditions A{Y,tj) — Aj{Y) for j = 1,2. After the gauge- 
transformation A' = A -I- V/, and using the above ex- 
pression for AS, we see that the probability amplitude 
connecting the transformed field configurations A'^ and 
A'2 is: 

^2'!' =^2iexp{z^(/(A2,/)-/(Ai,/))} (5) 

It is therefore natural to define the gauge-transformed 
wave-functional ^' by: 

*(^')-*(^)exp(i^/(A/)) (6) 

This definition ensures that '^{A,t) and '^{A,t) evolve 
according to the same set of probability amplitudes. 

In a Hamiltonian formulation, we associate to any clas- 
sical field configuration ^(r) a basis state |^). The 
gauge-transformation corresponding to / is now repre- 
sented by a unitary operator U (/) defined by: 

[/(/)|A)^cxp(*^/(A,/))|A') (7) 

The presence of the phase-factor is one of the essential 
features of the Chern-Simons term (i.e. the term propor- 
tional to v) added to the action. Note that when / varies, 
the family of operators U (/) gives rise to a representation 
of the full group of local gauge transformations. Indeed, 
at the classical level, the composition law in this group is 
given by the addition of the associated / functions, and 
because I{A, f) given in Eq. is itself gauge-invariant, 
we have U{f)U{g) = U{f + g). It is interesting to give 
an explicit expression for U{f). It reads: 

[/(/) = C/.=o(/) exp d^r (8) 

U,=o{f) = 1^ rf'r (SxH, + dyUy)/^ 

where 11 j; and 11^ are the canonically conjugated variables 
to Ax and Ay . Note that this no longer holds in the case 
of a manifold M with a boundary, as will be discussed in 
a few paragraphs. 

In the Hamiltonian quantization, a very important role 
is played by the gauge-invariant electric operators and 
Ey. In the absence of Chern-Simons term, they are sim- 
ply equal to and IIj,. When ^ 0, because of Eq. ((SJ, 
the transformation law for and IIj^ becomes: 

Ux ^Ux + vdyf 
IIj, By - vdxf 
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To compensate for this new gauge sensitivity of con- 
jugated momenta, the gauge-invariant electric field be- 
comes: 



Ey=ny + VA,, 

Any classical gauge-invariant Lagrangian gives rise, af- 
ter Legendre transformation, to a Hamiltonian which is 
a functional of E^, Ey, and B fields. If we add the 
Chern-Simons term to the original Lagrangian and per- 
form the Legendre transformation, we get a new Hamil- 
tonian which is expressed in terms of the new gauge- 
invariant E,j:, Ey, and B fields through the same func- 
tional as without the Chern-Simons term. For the special 
example of the Maxwell-Chern-Simons Lagrangian J^l, 
this functional reads: 

But although the Chern-Simons term preserves the 
Hamiltonian functional, it does modify the dynamical 
properties of the system through a modification of the 
basic commutation rules between E^ and Ey. More pre- 
cisely, we have: 

[E,{r),Eyir')]^-ihi2,.)S{r-r') (9) 



So it appears that finding the appropriate deformations 
of electric field operators plays a crucial role in construct- 
ing the Hamiltonian version of a Chern-Simons theory. 
We have also seen that such deformations are strongly 
constrained by the requirement of invariance under lo- 
cal gauge transformations. This discussion shows that 
most of the relevant information is implicitely encoded 
in the additional phase-factor I{A, f) involved in quan- 
tum gauge transformations, as in Eqs. © and 0. 

Let us now briefiy discuss what happens when the 
model is defined on a two-dimensional space manifold 
M with a boundary dM. Using Stokes' formula, we 
may recast the phase factor I{A, /) attached to a gauge- 
transformation as: 



I{AJ)^ [ d?vBf- [ /A.dl (10) 

Jm JdM 



In this situation, the phase factor I{A, /) is no longer 
gauge- invariant, and this implies that two gauge trans- 
formations attached to functions / and g do not commute 
because: 



iIiAJ) + IiA + Vf,g))-iIiA,g) + I{A + ygJ))= [ {fVg-gVf).dl (11) 

I 



In more mathematical terms, this reflects the fact that 
the phase-factor I{A,f), as used in Eqs. © and Q) 
defines only a projective representation of the classical 
gauge group, that is: 

U{f)U{g) = cxp (^-z^ /Vg.dl^ U{f + g) (12) 

As first shown by Witten2&, this may be understood in 
terms of a chiral matter field attached to the boundary of 
M . An explicit example of boundary degrees of freedom 
induced by a Chern-Simons term has been discussed re- 
cently in the case of a Z2 model on a triangular latticaSl. 

To close this preliminary section, it is useful to discuss 
the case of a finite cylic group 'S.n- In the U{1) case, for 
a pair of points r and r', we have a natural group element 

defined by exp(i-^ A.dl) where the integral is taken 
along the segment joining r and r', and $0 is the flux 
quantum in the model. For a finite group G, the notion 
of a Lie algebra is no longer available, so it is natural to 
define the model on a lattice. In a classical gauge theory, 
each oriented link ij carries a group element gij G G. We 



have the important constraint gijgji — e, where e is the 
neutral element of the group G. In the quantum version, 
the Hilbert space attached to link ij is the finite dimen- 
sional space generated by the orthogonal basis \gij) where 
gij runs over all elements of G. For a lattice, the corre- 
sponding Hilbert space is obtained by taking the tensor 
product of all these finite dimensional spaces associated 
to links. In the model, gij becomes an integer mod- 
ulo N, pij. The connection with the continuous case is 

obtained through the identification J^"* A.dl = j^Pij- On 
each link ij, we introduce the unitary operator tt^ which 
sends \pij) into \pij + 1). In the absence of a Chern- 
Simons term, the generator of the gauge transformation 
based at site i (which turns pjk into pjk + Sji — Ski) is 

Ui — Hj'^ ""iji where the product involves all the nearest 
neighbors of site i. By analogy with the continuous case, 
the presence of a Chern-Simons term is manifested by an 
additional phase-factor whose precise value depends on 
the lattice geometry and is to some extent arbitrary, since 
fluxes are defined on plaquettes, not on lattice sites. On 
a square lattice, a natural choice is to define Ui according 
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09. 



(i) 



J2 pj^) 



(13) 



where C{i) is the oriented loop defined by the outer 
boundary of the four elementary plaquettes adjacent to 
site i. This expression has exactly the same structure as 
Eq. but somehow, the local magnetic field at site i is 
replaced by a smeared field on the smallest available loop 
centered around i. It has been showrkiS that a consistent 
quantum theory can be constructed only when ly/h is an 
integer multiple of N/tt. 



III. GENERATORS OF LOCAL GAUGE 
TRANSFORMATIONS 

As discussed in the previous section, the most impor- 
tant step is to construct a non-trivial phase factor which 
appears in the definition of unitary operators associated 
to local gauge transformations, generalizing Eq. Q. For 
this, let us first define the operator Lij{g) which is the left 
multiplication of gij by g, namely: Lij{g)\gij) = \ggij). 
For any site i and group element g, we choose the gener- 
ator of a local gauge transformation based at i to be of 
the following form: 

U^i9)^llL,,{g)YlxH^M9) (14) 

where j denotes any nearest neighbor of i and $(i,r) is 
the flux around any of the four square plaquettes, cen- 
tered at r, adjacent to i. Here, and in all this paper, 
we shall focus on the square lattice geometry, to simplify 
the presentation. But adaptations of the basic construc- 
tion to other lattices are clearly possible. Since we are 
dealing with a non-Abelian group, we have to specify an 
origin in order to define these fluxes, and it is natural to 
choose the same site i, which is expressed through the 
notation $(i,r). Since we wish Ui{g) to be unitary, we 
require |x$(5)| = 1- It is clear from this construction 
that two generators Ui{g) and Uj{h) based on different 
sites commute, since the phase factors x*(5) ^-re gauge 
invariant. 

This form is a simple generalization of the lattice 
Chern-Simons models for the cyclic groups 'Kn discussed 
in the previous section. In this example, for a square pla- 
quette ijkl, the flux <& is equal to pij +pjk +Pki +Pii mod- 
ulo N, and g is simply any integer modulo N. Eq. H13|l 
above may be interpreted as: 



X$(g) = exp 



(15) 



This is a well defined function for $ and g modulo N 
only if ly/h is an integer multiple of 2N /-k. We have not 
succeeded to cast odd integer multiples of N/tt for v/hva 



the framework of the general construction to be presented 
now. This is not too surprising since these models were 
obtained by imposing special periodicity conditions on 
an infinite-dimensional Hilbert space where pij can take 
any integer valuc^^. Our goal here is not to write down 
all possible Chern-Simons theories with a finite group, 
but to easily construct a large number of them, therefore 
allowing for non-trivial phase-factors when two localized 
flux excitations are exchanged. 

As discussed in the Introduction, a very desirable prop- 
erty, at least for the sake of finding possible physical im- 
plementations, is that these deformed generators define 
a unitary representation of the group G. So we wish to 
impose: 



U^{g)U,{h) = U,{gh) 



or equivalently: 



X/i*/i-i(5)x*(^) = X'i'igh) 



(16) 



(17) 



To solve these equations let us first choose a group el- 
ement Let us denote by the stabilizor of <f> un- 
der the operation of conjugacy, namely h belongs to 
whenever h^h~^ = <&. This notion is useful to describe 
the elements of the conjugacy class of $. Indeed, we note 
that gh^{gh)^^ = g^g~^ if h belongs to Hg,. Therefore, 
the elements in the conjugacy class of $ are in one to 
one correspondence with the left cosets of the form gH^. 
Let us pick one representative in each of these cosets. 
We shall now find all the functions Xg ^g-^ig)- First we 
may specialize Eq. (|17(1 to the case where h belongs to 
H^, giving: 



X*(5)X<I>(^) = X<i'{gh) 



(18) 



In particular, it shows that the function h X$(^) de- 
fines a group homomorphism from to U{1). Once 
this homomorphism is known, we can specify completely 
X<s>ig) for any group element g once the values X*(5n) for 
the coset representatives are known. More explicitely, we 
have: 



Xi-ignh) = X'i>(5n)x*(^) 
where h e Finally, Eq. ltT7|l yields: 

, N _ Xf(gff«) 



(19) 



(20) 



Let us now show that for any choice of homomorphism 
h X*(^)7 ^ S and unit complex numbers for 

X*(gn), Eqs. ((T^ . (Pn|l reconstruct a function x*(3) 
which satisfies the condition H17|) . Any element g' in G 
may be written as g' — gnh, with h G We have: 



Xg"s>g'-^{g) = Xg^^g-^{g) 



m X*(35rO IB 



X'i'ign) 

x$(gg«fa) ^ X'fjgg') 
x^ignh) Xfig') 



(21) 
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which is exactly Eq. (|17|l . 

Note that there are many equivalent ways to choose 
these functions x*(5)- Let us multiply the system wave- 
function by a phase-factor of the form JJ- Jj'^^ e(^(«, r)), 
where |e(^)| = 1- Under this unitary transformation, 
X*(5) is changed into xi'ig) given by: 



(22) 



In particular, it is possible to choose the values of 
e((ji„$g^^) so that x<s>{gn) = 1- Although this does not 
seem to be required at this stage of the construction, it 
is also necessary to assume that when the flux $ is equal 
to the neutral element e, Xe{g) = 1- This will play an 
important role later in ensuring that the phase-factor ac- 
cumulated by the system wave-function as a fluxon winds 
around another is well defined. 



IV. 



BASIC PROCESSES FOR FLUXON 
DYNAMICS 

A. General description 
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FIG. 1; The site labeling convention. For any bond (ij) 
(shown here as middle vertical bond) the surrounding sites are 
labeled 1, 2, 3, 4 as indicated here. The fluxes $(i, L), $(i, R), 
{^{j, L), ^{j, R)) are counted counterclockwise, starting from 
site i, (j), e.g. $(i,L) = grjgjmidu 

Our goal here is to construct local gauge invariant op- 
erations for basic fluxon processes: fluxon motion, cre- 
ation of a pair with vacuum quantum numbers, branch- 
ing of a single fluxon into two fluxons and their time 
reversions. These three types of elementary processes 
can all be derived from a single operation, the electric 
field operator that at the level of classical configura- 
tions, is simply a left multiplication Lij{g) attached to 
any link ij. To show this, let us consider a pair of two 
adjacent plaquettes as shown on Fig. ^ We denote by 
<i?(«,L) (resp. <l>(i,i?)) the local flux through the left 
(resp. right) plaquette, with site i chosen as origin. Sim- 
ilarly, we define fluxes ^{j,L) and ^{j,R). Changing 
the origin from i to j simply conjugates fluxes, according 
to $(j, L(i?)) = gji^{i, L{R))gij. The left multiplication 
Lij{g) changes gij into g'^j = ggij. Therefore, it changes 



simultaneously both fluxes on the left and right plaque- 
ttes adjacent to link ij. More specifically, we have: 



^'{i,L)^g^i,L) 
^'ii,R)^^i,R)g-^ 



In particular, this implies: 



$'(i, R)^'ii, L) = $(i, i?)$(i, L) 
^'{i,L)^'{i,R)^g^{i,L)^{i,R)g~^ 



(23) 
(24) 

(25) 
(26) 



Note that transformation laws for fluxes based at site j 
are slightly more complicated since they read: 



$'(j,L) = H3,L){gj,gg,j) 
^'{j,R)^{gy,g~^g,,)Hj,R) 



(27) 
(28) 



This asymmetry between i and j arises because 
g'ji = gjig~^^ so we have: 



Lijig) = R]i{g ^) 



(29) 



where Rji{h) denotes the right multiplication of gji by 
the group element h. In the absence of Chern-Simons 
term, Lij{g) commutes with all local gauge generators 
with the exception of Ui{h) since: 



U,{h)L,,{g)Udh)-^ - L,j{hgh-^) 



(30) 



We now apply these general formulas to the elementary 
processes involving fluxes. Suppose that initially a flux 
<i> was localized on the left plaquette, and that the right 
plaquette is fluxless. Applying Lij{g — <i>~^) on such 
configuration gives: 



$'(i,L) = e 
$'(i,i?) $ 



(31) 
(32) 



which shows that a $-fluxon has moved from the left to 
the right plaquette. A second interesting situation occurs 
when both plaquettes are initially fluxless. The action 
of L.y((f>^^) on such state produces a new configuration 
characterized by: 



<i>'{i,L) ^ <i>-^ 
$'(i,i?) = $ 



(33) 
(34) 



So we have simply created a fluxon and antifluxon pair 
from the vacuum. Of course, applying Ly($) on the 
final state annihilates this pair. Finally, a single flux 
$ = <I>2$i originaly located on the left plaquette may 
split into a pair $i on the left and $2 on the right. This 
is achieved simply by applying Lij{^2^)- 

In order to incorporate these elementary processes into 
a Hamiltonian Chern-Simons theory, we have to modify 
Lij {g) into an electric field operator f ^ (g) by introducing 
phase-factors so that it commutes for all generators Uk{h) 
ii k ^ i and that: 



U,ih)£,jig)U,ih)~' = S^jihgh- 



(35) 
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As explained in the introduction, we shall need 
only to construct £ij{g) for special types of config- 
urations for which at least one of the four fluxes 
L), ^(i, i?), i), <j>'(i, i?) vanishes. Nevertheless, 
it is useful to first construct Eij (g) for an arbitrary flux 
background. The requirement of local gauge-symmetry 
induces strong constraints on the phase- factors as 
we shall now show. These constraints are less stringent 
when we restrict the action of £ij (g) to the limited set of 
configurations just discussed. 



Construction of gauge-invariant electric field 
operators 




hj h 
Generic initial state 

FIG. 2: Construction of the gauge-invariant electric field op- 
erator £ij(g). This operator transforms the pair of fluxes 
($L, ^r) of the left and right plaquettes into the pair 
((?$L, ^Rg~^). We define the amplitude of the transition 
between two reference states shown at the top of the fig- 
ure by j4($l, g). The amplitude of the process start- 
ing from a generic initial state shown in the lower left 
can be related to A{^L,^R,g) using the gauge invariance. 
This is done in two steps. First, the gauge transformation 
Ui{h}) = U^_{hi)U2{h2)Ui{hs)UA[hA)Uj{hj) is used to relate 
the amplitude of the process starting from the generic state to 
the amplitude of the transition a(<&L, ^R,g) between reference 
state and a special state shown in the middle right. Second, 
we use gauge transformation W{g) = U4,{g)Uj{g)Uz{g) to re- 
late the special state to the reference state on the upper right. 
Site labeling is the same as in Figure 1. 

We now construct the electric field operators £ij {g) at- 
tached to links. In the absence of Chern-Simons term the 
electric field operator is equivalent to the left multiplica- 
tion of the link variable: 



(36) 



In this case, £ij{g) commutes with all local gauge genera- 



tors with the exception of Ui{h), and Eq. (|35|l is satisfied. 
We have also the group property, namely: 



(37) 



The Chern-Simons term gives phase factors x*(5) fo 
the gauge generators, so we expect some phase factor, 
Tij, to appear in the electric field operators as well: 



£ 



Li'j Tj 



The gauge invariance condition allows us 



to relate the phase factors associated with different field 
configurations to each other. Specifically, we introduce 
the reference states (shown on the top of Fig. |2 ) in 
which only two bonds carry non-trivial group elements. 
We define the transition amplitude induced by the elec- 
tric field between these reference states by A{^l, g). 
In order to find the phases Tij for arbitrary field config- 
uration we first transform both the initial and final state 
by U{{h}) = Ui{hi)U2ih2)U3ih3)Ui{hi)Uj{h,). This re- 
lates the amplitude of the generic process to the ampli- 
tude, a{^L,^R, g), of the process that starts with the 
reference state but leads to the special final state shown 
in Fig. El middle right. Collecting the phase factors asso- 
ciated with the gauge transformation U{{h}) before and 
after the electric field moved the flux we get 

^ _ X<s>'ii,L){hi) X^'(i,R){h2) X<i'{j,R){h:i) 

X<[.(i,L)(^l) X<I.(i,fi)(/l2) X^u,R){hz) 
^ X$'(j,L)(^4) X$'(j,L)(^j) Xj"(j\i?.)(^j) 
X*(j,L)(^4) X*(j,L)(^j) X'i>(]M){hj) 



X 



a($L,$i?,.g) 



(38) 



where $(«,L) denotes the flux in the left plaquette 
counted from site $(j, R) denotes flux in the right 
plaquette counted from site j, and prime refers to the 
flux configuration after the action of the electric field. 
Finally, we employ the gauge transformation W{g) — 
Ui{g)Uj{g)U^{g) to relate this special final state to the 
reference state. The phase factor associated with this 

gauge transformation is (x*i,g(5)Xg-i*ji (5)) so 



a{<^L,^R,g) 



A{<^L,^R,g) 



(x*i.s(5)Xs-i<i.«(.9))' 



(39) 



In order to express the phase factors, Tij, through the 
initial field configuration we relate the parameters, /ifc, of 
the gauge transformations to the bond variables gk by 



hi = gu, h2 = g2i, 
hi — gAjgji, hj ~ gji 



93jgji, 



and the fluxes in the left and right plaquettes before and 
after electric field operator has changed them. Before the 
electric field operator has acted the fluxes were 

*(*,i) = $L 

^{j,L) = gJi^Lg^j 

$(l,i?) = 

^(j, R) = 93i^R9ij 
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while afterwards they become 

Combining the preceding equations and using the rela- 
tion H17() a few times, we get the final expression for the 
phase-factor, 

^ _ Xl>' {i,L){gii) X<S>'{i,R}{g2i) X<S>'U,L){94j) 
X*(i,L)(5li) X$(i,fl)(fl'2i) X$(i,L)(ff4j) 

^ X$'(j\fl)(.93j) X<s>'(i,L){9ji)'^ X<s>'{i,R){9'jtY 

X<f(j,R)i93j) X<S>{i,L){9]i)^ X'S>(i,R){9ji? 

x^($i,$H,g) (40) 

where we have used the definition g[j — ggij 
{9'ji — 9ji9~^) to make it more symmetric. 

Commutation of £ij{g) with C/2(/i2), U3{h3), 

J74(/i4), and Uj{hj) follows directly from this construc- 
tion. It can also be checked directly from (|4()|l . using 
the condition H17|l on the elementary phase- factors x*(.9)- 
Note that sites i and j play different roles, which is ex- 
pected because Sijig) acts by left multiplication on gij 
whereas Sjiig) acts by right multiplication on the same 
quantity. 

Although the electric field operator commutes with 
Ui{hi), U2{h2), (/^{hs), Ui{hi), and Uj{hj), it does not 
necessarily commute with Ui(h) even if hgh~^ — g. In 
fact, the requirement of this commutation leads to an im- 
portant constraint on possible choices of phases X*(5)- 
The appearance of the new constraints becomes clear 
when one considers a special field configuration shown 
in Fig. 121 Two identical field configurations shown on 
the top and bottom left of this figure can be obtained by 
two different gauge transformations from the reference 
state if both $i and commute with h: in one case 
one applies gauge transformation only at site i , in the 
other one applies gauge transformation on all sites ex- 
cept i. Provided that the resulting states are the same, 
i.e. gh — hg, the total phase factor Tij obtained by these 
two different ways should be the same. 

The phase factors associated with these gauge- 
transformations are the following 

Ui{h)-^Xi>L{h)x-fR{h) 

Ur\h)^ TTT^ TTT 

XgfL(h)X<S>ng-^[h) 

^'^^^^xlJh-^)Xg^Ah-')xl^^Jh-^)x^^,-.ih-^) 
Putting all these factors together we conclude that the 



fiii(g) 




Ui(h) 



U:{h-1) 



6ii{g) 



FIG. 3: (Color online) Gauge invariance of electric field oper- 
ator implies condition on the phase factors X*(5)- The state 
shown at the top and botton left of the figure can be trans- 
formed to the reference state (middle) in two ways: by apply- 
ing a gauge transformation V — Ylky^i^i^W shaded 
sites or by applying transformation Ui{h^^) on site i. If 
c&Lft = ti^L and — fi^R these transformations lead 

to the same reference state shown in the middle. Further- 
more, the action of electric field on these states can be found 
by making a gauge transformation of the final state shown 
on the middle right if gh = hg. The phases of electric field 
operator obtained in these two ways should be equal giving 
additional condition on the phase associated with the gauge 
field transformation. 



gauge invariance of the electric field operator implies that 

X$i, Wx$fl (fe) ^ 
Xg*i(/i)x*flff-i(M 

xlAh-^)xlJh-^) 

xlr,gih~^)Xg'fL{h-^)xl-i^Jh-^)X^^g-i{h-^) 



(41) 



This condition can be further simplified by using the 
main phase factor equation H17|) . We start by not- 
ing that because h commutes with ^r, and g, 

X<s>H,Lih)x<fR.L{h~'^) = 1 and Xs-Ch,!, (MXai-^.i (^"^) = 1- 
This gives 

xlr,i^)xln(^) = xl^g{h)xl<s>^ih)xl-i<s>nih)xlj,g-i{h) 



Furthermore, combining the identities 

Xg<fLih) = Xg*I.(5/^.9"^) 

= Xi>Lgi9)x<S'Lg{h)Xg<S'Li9~^) 

and 

1 = Xg^Ligg^^) = X<S'Lgi9)Xg<S>L{9~^) 



(42) 
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we get 



(43) 
(44) 



This reduces the condition (|42|l to a much simpler final 
form 



X'S'lW X$r(^) 



= 1 



(45) 



We emphasize that constraint l|45|l on the phase factors 
has to be satisfied only for fluxes satisfying the condi- 
tion {h^Lh~^,h^Rh~^,hgh~^) = ($L,$fl, g). Although 
we have derived this condition imposing only the gauge 
invariance of the electrical field acting on a very special 
field configuration, a more lengthy analysis shows that it 
is sufficient to ensure that in a general case 



U,{h%,ig)\^) = £,,ihgh-')U,{hm 



(46) 



The details of the proof are presented in Appendix ^ 

Unlike Eq. (|17|l . the constraint (|45(l relates functions 
X$(g) and x*'(5) for ^ ^^id $' belonging to different 
conjugacy classes. As shown in section IVI Bl this con- 
straint strongly reduces the number of possible Chern- 
Simons theories. Note that if both ^^^''(g) and xi^\g) 
arc solutions of the fundamental relations (|17l) and (|45|l . 

their product xi^'' (5)x$ (s) is also a solution. So there 
is a natural group structure on the set of possible Chern- 
Simons models based on the group G, which is transpar- 
ent in the path-integral description: this means that the 
sum of two Chern-Simons action is also a valid Chern- 
Simons action. 

Is this construction also compatible with the group 
property 137|l ? From Eq. (|45l) . we obtain: 



where: 



/3($L,$fl,5) 



A{^L:^R,9)A{g^L,^B,g~ 
A{<^L,^R,g'g) 



It does not seem that (3{^l,^r, g) are always equal to 
unity for any choice of x* {g) that in turns determines the 
amplitudes A{s) (see Appendix IXjl . But this is not re- 
ally a problem because this group property does not play 
much role in the construction of gauge-invariant Hamil- 
tonians. 

We now specialize the most general constraint arising 
from local gauge-invariance to the various physical pro- 
cesses which are required for fluxon dynamics. For the 
single fluxon moving operation, we have $l = <I>, = e, 
and g = so the condition Eq. (|45|l is always satisfled. 
For the pair creation process, we have = ^r — e, and 
g = The constraint becomes: 



Finally, let us consider the splitting of an isolated fluxon 
into two nearby ones. This is described by = ^2^1, 
<f>_R = e, and g — ^2^- We need then to impose: 



Xi>i{h)X'fAh )Y 



(x*-i(^)x*(^)) 



1 if h<Ph- 



(48) 



= 1 if h^2h^'^) = ($1, $2) 

(49) 

It is clear that condition H48|) is a special case of 
the stronger condition (|49|l . Furthermore, multiplying 
the conditions H49|l for pairs of fluxes {^l,^r) and 
(.g$L, $7?(7^^) we get the most general condition H45|l : 
this shows that constraint (|49|) is necessary and suffl- 
cient condition for the gauge invariant deflnition of elec- 
tric fleld operator acting on any flux configuration. 



V. GENERAL EXPRESSION FOR HOLONOMY 
OF FLUXONS 



Let us consider two isolated fiuxcs carrying group ele- 
ments gi and g2, and move the first one counterclockwise 
around the second one, as shown on Fig. 0| This can be 
done by successive applications of local gauge-invariant 
electric field operators as discussed in the previous sec- 
tion. Although we wish to work in the gauge-invariant 
subspace, it is very convenient to use special configura- 
tions of link variables to illustrate what happens in such 
a process. We simply have to project all special states 
on the gauge invariant subspace, which is straightforward 
since the fluxon moving operator commutes with this pro- 
jector. The initial fluxes are described by two vertical 
strings of links carrying the group elements gi and g2, 
see Fig. 01 When several other fluxes are present, besides 
the two to be exchanged, it is necessary to choose the 
location of the strings in such a way that no other flux 
is present in the vertical band delimited by the strings 
attached to the two fluxons. During the exchange pro- 
cess, the first fluxon collides with the second string, and 
after this event, it no longer carries the group element gi, 
but its conjugate g[ = g2^gig2- After the process is com- 
pleted, the flrst flux has recovered its original position, 
but the conflguration of group elements has changed. If 
we measure the second flux using a path starting at point 
O shown on Fig.0|we find that it has also changed into 
g'2 — 52" ^5r^ 5231 32- The final state can be reduced to its 
template state built from two strings carrying groups el- 
ements g'l and g'2 by the gauge transformation Y[i Ui{hi) 
where hi is locally constant in the two following regions: 
the core region inside the circuit followed by the first 
fluxon {hi = /icore), the intermediate region delimited by 
the two initial vertical strings and the upper part of the 
circuit [hi = /tint)- Note that because we do not wish to 
modify external fluxes, we cannot perform gauge trans- 
formations in the bulk, outside of these regions. 

Group elements /icore and /imt have to satisfy the fol- 
lowing conditions, which may be obtained readily upon 
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VI. APPLICATION TO DIHEDRAL GROUPS 







if 








—I- 


uu: 






c 





FIG. 4: (Color online) Braiding of two fluxes yields in a 
non-trivial transformation of their values and a phase fac- 
tor. We start with the flux configuration shown in upper 
pane with fluxes gi (right) and 172 (left) connected to the 
edge by the strings of gi (purple arrow) or g2 (black) group 
elements. Moving the right flux around the left leaves be- 
hind a string of gi elements until its path crosses the vertical 
string of g2- Upon crossing the string changes to the string 
oi g'l = 92^9192 (rsd arrows). Performing the gauge transfor- 
mations with h — (g'l) on all sites indicated by full (red) dots 
and with h = gih on sites indicated by empty (orange) dots 
reduces the configuration to the template shown in the last 
pane with the new fluxes g'l and g'2 (see text, Eq. 1521 '). 



inspection of Fig. 0| 



^corc 


/ 

= .9i 


7 /-I 
'^corc^l 


= e 




= .92 


h 

'^corc"'int 


= 31 






^corc^2^corc 


= 92 



These equations are mutually compatible, and we get 
^int = 91^92^9192- Since fluxes are present in the core 
region, they will contribute a phase- factor /(51, 52) when 
the gauge transformation from the template to the actual 
final state is performed. The final result may be summa- 
rized as follows: 



lffl,52) 

.9'i 



52)15^,52) 



.92 ^9192 



92 = 92 ^9i ^929192 
' ^2, 



/(5i,52) = Xg'Sa'i) Xs^(.9'i) 



(50) 
(51) 
(52) 
(53) 



A. General properties of dihedral groups 

The dihedral groups are among the most natural 
to consider, since they contain a normal cyclic group 2jv 
which is simply extended by a ^2 factor to form a semi- 
direct product. In this respect, one may view this family 
as the most weakly non-Abelian groups. can be de- 
scribed as the isometry group of a planar regular polygon 
with TV vertices. The Zjv subgroup corresponds to rota- 
tions with angles multiples of 2n/N . We shall denote by 
C the generator of this subgroup, so C may be identified 
with the 2ti/N rotation. Dm contains also N reflections, 
of the form rC". The two elements C and r generate 
_Djv, and they are subjected to the following minimal set 
of relations: 



tCt = C- 



(54) 
(55) 
(56) 



This last relation shows that indeed Dn is non-Abelian. 

The next useful information about these groups is the 
list of conjugacy classes. If = 2M + 1, contains 
M + 2 classes which are: {e}, {C,C-i},..., {C*^C-^^}, 
{t,tC, ...,tC^-^}. If N = 2M, there are M + 3 classes 
whose list is: {e}, {C^^}, {C,C-^},..., {C^^-\C-^^+^}, 
{r,rC^...,rC^-'},{rC,TC^...,TC^-'}. 

As shown in section IIIII in order to construct possi- 
ble phase factors associated to gauge transformations, 
we need to know the stabilizors of group elements for the 
conjugacy operation of Djv acting on itself. Here is a list 
of these stabilizors: 
For N odd: 



e 
CP 

TCP 



Dn 

1 < p < - 1 
{e,rCP}, 0<p<N 



For N even: 

e 

CN/2 



CP 
tCP 



Dn 



1<P<N-1, pj^ 



N 



{e,C^/2,TCf,rCP+^/2|^ 0<p<Af-l 



The new phase-factor /(f/i, (72) appears not to depend on 
the detailed path taken by the first fluxon, but just on the 
fact that it winds exactly once around the second. In this 
sense, our construction really implements a topological 
theory. 



Finally, we need to choose homomorphisms from these 
stabilizors into U{1). In the case of a cyclic group '2.N 
generated by C, homomorphisms x are completely de- 
termined by x(C), so that x(C^) = xi^Y: "^ith the con- 
straint: x(C)^ = 1. For the group Dn itself, we have: 
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x{t'^CP) = xi'TYxi^y^ with the following constraints: 

X(C)^ = 1 (57) 
X(r)2 = 1 (58) 



1 



(59) 



These are direct consequences of generator rela- 
tions (|5Hl . (|55|l . l(K^ . Again, the parity of N is relevant. 
For TV odd, x(C) — 1, which leaves only two homomor- 
phisms from into U{1). For N even, x{C) = ±1, 
and there are four such homomorphisms. The last pos- 
sible stabilizor to consider is the four element subgroup 
of D2M, S = {e,C^\TCP,TCP+^'}. This abelian group 
has four possible homomorphisms into U{1), which are 
characterized as follows: 

X(C'0 = ±1 (60) 
xirCn - ±1 (61) 

x{rCP+^') = xirCnxiC'') (62) 

B. Classification of possible models 

1. N odd 

Let us first consider conjugacy classes of the form 
{CP, C^P}. Since the stabilizor of for the conjugacy ac- 
tion of Dn is Zat, we have: xcp{C'^) = ^p, where ujp = 1. 
Choosing Xcp(t) = 1, we have for fluxes in the cyclic 
group generated by C: 



xc4n = Lo^p 



N 



UJ, 



UJpUJ—p — 1 



(63) 
(64) 
(65) 
(66) 



For the remaining conjugacy class, we have XriT) = rj, 
with J] — ±1. Choosing Xr{CP) = 1, and using Eqs. H19() 
and l(2Il, we obtain: 



xrc^n^i 



(67) 
(68) 



All these possible phase-factors satisfy the following 
property: 



X<i.-i(ft-)x*(^) = 1 



(69) 



so that Eq. 148|) is always satisfied. So no new constraint 
is imposed by the requirement to create or annihilate 
a pair of fluxons. What about the stronger condition 
Eq. (|49|l ? Its form suggests that we should first deter- 
mine pairs of fluxes ($i,$2) such that their stabilizors 
-ff$j^ and H^^ have a non-trivial intersection. This occurs 
if both <f>i and $2 are in the Z^v normal subgroup gener- 
ated by C, or if $1 = $2 = tC^. The second case simply 
implies Xt{t)'^ = 1, which is not a new condition. In the 



first case, choosing h in Za? as well shows that x*(^)^ 
is a homomorphism from Zjy to U(l) with respect to 
both $ and h. This is satisfied in particular if x$(^) it- 
self is a group homomorphism in both arguments. This 
sufficient (but possibly not necessary) condition simpli- 
fies algebraic considerations; it can be also justified from 
physical argument that the theory should allow for the 
sites with different number of neighbours, Z, which would 
change x-fW'^ ^ X<f{hY- 

Replacing this the constraint on x*(^)^ by the con- 
straint on x<s>{h)^ we get 



XcAn^ujP 



c.^-1 



(70) 
(71) 
(72) 



Therefore, the class of possible phase-factors (which is 
stable under multiplication) is isomorphic to the group 
Zat X Z2. This group of phase- factors is identical to 
the group of possible Chern-Simons actions since in this 
case, H'^{Dn,U{1)) = Zat x Z^Si. Very likely, the co- 
incidence of these two results can be traced to the ab- 
sence of projective representations for D^r for N odd^-: 
as explained above, the projective representations are not 
allowed in our construction but are allowed in the classi- 
fication of all possible Chern-Simons actionsSi. 



2. N even 

The conjugacy classes of the form {C^,C^p} behave in 
the same way as for N odd. So writing N — 2AI, we 
have: 



Xcp(C'') = c.^ 



UJ, 



N 



p i {0,M},mod 2iV 



(73) 
(74) 
(75) 
(76) 
(77) 



The conjugacy class {C*^} is special since its stabilizor is 
Z?Ar itself. As discussed in Section IVl Al above, there are 
four homomorphisms form Dj^ into C/(l) that we denote 
by: 



uj,ujm e {1, -1} 



(78) 
(79) 



Let us now turn to Xrid) with the corresponding stabi- 
lizor equal to {e,C^^ ,t,tC^^}. As seen in section IVl Al 
the four possible homomorphisms may be written as: 



Xr(T) = ryo e {±1} 

XriC^') = Co e {±1} 



(80) 
(81) 
(82) 
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From this, we derive the expression of Xrig), in the fol- 
lowing form (0 < p < M — 1): 



XrirC-n = Vo 



(83) 
(84) 
(85) 
(86) 



Furthermore, because C^rC p = tC ^p, Eq. 1201) imphes: 



Xrc-^Ag) = XrigCP) 



(87) 



The last conjugacy class to consider contains rC, with 
the stabihzor {e,C^^ ,tC,tC^^^^}. In this case, we may 
set (0 < p < M - 1): 



Xrcin = 1 

Xrc(C^+'0 = Ci e {±1} 
XrcirC'-p) = m e {±1} 

Xrc^-^Ag) = XrcigCP) 



(89) 
(90) 
(91) 
(92) 



Here again, the constraint Eq. H48|l is always satisfied. 
To impose Eq. H49|l . we have to consider pairs of fluxes 
(<i>i,$2) such that n -ff$2 is non-trivial. As before, 
choosing $i and <I>2 in the Zjv subgroup imposes tUp = 
ui'P, with uj^ = 1. A new constraint arises by choosing 
$1 = CP and $2 = tCp' . In this case, C^^ belongs to 
their common stabilizor. Eq. (|49|) implies: 

XcAC'')XrcAC'') = Xc.rcAC') = XrC'^AC'') 

But using Eq. 117() this yields: 



(93) 



Therefore we have the constraint: Co = Ci = 1; which en- 
ables us to simplify drastically the above expression for 
phase-factors: 



and 



XrcAC) = 1 
XrC^piTC^) = VO 



X.C^P+i(C«) = l 
XrC2p+i(TC'') = T]i 



(95) 

(96) 



(97) 
(98) 



But Eq. H93|l now implies that Xcp(C*^) = 1 for any p, 
which is satisfied only when = 1. Specializing to 
p — M , we see that the common stabilizor contains now 
two more elements, namely rC^ and rC+^. Eq. H93(l 
now requires that: 



M even 
Modd 



Co ^ 1 



(99) 
(100) 



It is then easy to check that considering the common sta- 
bilizor of tCp and rC , which may contain two or four el- 
ements, does not bring any new constraint. Finally, w be- 
longs to the group and among the three binary vari- 
ables rjQ, rji, and uj, only two arc independent. Therefore, 
the set all all possible phase-factors for the D2A1 group 
is identical to Z7\/ x '2.2 x '2.2. This contains only half 
of H^{D2M,U{1)) which is equal to Z2M x Z2 x Zi^"*. 
But D2M admits non trivial projective representations 
since H"^ {D2M ,U {!)) = Z2 which cannot appear in our 
construction. The important result is that in spite of 
this restriction, we get a non-trivial subset of the possi- 
ble theories also in this case. 



C. Holonomy properties 

Using the general expression 1)53(1 . and the above de- 
scription of possible phase-factors, we may compute the 
adiabatic phase induced by a process where a fluxon gi 
winds once around another fluxon 32- The results are 
listed in the last column of table HI and they are valid for 
both parities of N . 



TABLE I: Adiabatic phase f{gi,g2) generated in the pro- 
cess where fluxon gi winds once around a fluxon g2. Values 
of phase- factors Xg[{g'i) ^-nd Xg'^ig'i) ^re given for dihedral 
groups Dn with odd A^. For even N, expressions for Xg[ (Si) 
and Xg' (gi) s-re slightly more complicated, but interestingly, 
the main result for /(<?i, (72) is the same as for odd N, namely 
it involves only the complex number u. 
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VII. CONCLUSION 



Generally, in order to appear as a low energy sector 
of some physical Hamiltonian, the Chcrn-Simons gauge 
theory has to involve gauge transformations that depend 
only on a local flux configuration. Furthermore, to be in- 
teresting from the view point of quantum computation, 
the theory should allow for a local gauge invariant elec- 
tric field operator that moves a flux or fuses two fluxes to- 
gether. Here we have analyzed non-Abelian gauge mod- 
els that satisfy these general conditions; our main result 
is the equation for the phase factor x associated 
with a gauge transformation. Furthermore, we have com- 
puted the flux braiding properties for a given phase factor 
that satisfies these conditions. Finally, we have applied 
our general results to the simplest class of non-Abelian 



13 



groups, dihedral groups D„. The fluxon braiding prop- 
erties in these groups are summarized in Table I. 

Inspection of the Table I shows that even for the small- 
est groups, Chern-Simons term modifies braiding proper- 
ties in a very non-trivial manner and gives phases to the 
braiding that were trivial in the absence of the Chern- 
Simons term. In the scheme^! where the pair of two 
compensating fluxes {tCp,tCp) are used to encode one 
bit, the transformations allowed in the absence of Chern- 
Simons term are limited to conjugation. In the pres- 
ence of Chern-Simons term the braiding of such bit with 
the controlled flux results in a richer set of transforma- 
tions that involve both conjugation by group elements 
and phase factors (see Table I) but does not change the 
state of the flux as it should. We hope that this will make 
it possible to construct the universal quantum computa- 
tion with the simplest group that does not involve 
operations that are difficult to protect (such as charge 
doublets). 

The implementation of the microscopic Hamiltonians 
discussed in this paper in a realistic physical system is 
a challenging problem. The straightforward implemen- 
tation would mean that the dominant term in the mi- 
croscopic Hamiltonian is H = —t^igUi{g) so that all 
low energy states are gauge invariant. This is not easy 
to realize in a physical system because operator U{g) in- 
volves a significant number of surrounding bonds. We 
hope, however, this can be achieved by a mapping to 
the appropriate spin model as was the case for Abelian 
Chern-Simons theories; this is the subject of the future 
research. 
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FIG. 5: (Color online) Derivation of the general condition on 
the phase factors from the requirement of electric gauge field 
invariance. The action of the electric field on arbitrary states 
(middle of the figure) related to each other by the gauge trans- 
formation at site i (shown as black dot) can be expressed via 
the phase factors A{s) and A{h{s)) characterizing the transi- 
tion between reference states (bottom and top). This involves 
gauge transformations on all grey (blue) sites. The gauge in- 
variance implies that the total phase factor accumulated when 
moving around the diagram is unity which implies the con- 
dition on the ratio A{s) / A{h{s)) . As in text, s stands for a 
triple ($i,$i{,(;), and h{s) is obtained by conjugating the 
components of s by h. 



APPENDIX A: ELECTRIC FIELD GAUGE 
INVARIANCE: GENERAL CRITERION 

In the absence of additional constraints on the gauge 
transformation phase factors, X$(5), the definition of the 
electric field operator implies that it is gauge invariant 
under gauge transformations on all sites except i while 
the the gauge transformation on the latter gives 

U,ih)£,,{g)\'<S/) ^ ai^L,^R,g\h)£^,ihgh-')Udh)\^) 

(Al) 

(A2) 
(A3) 



with 



a(*L,$i?,5lM 



A{<pL,^B.,g) 

A{h<^>Lh-\h^Rh-\hgh-^) 



Here the notations are those of Fig. 13 where state 
corresponds to the initial state on the figure. The compu- 
tation leading to this equation is depicted on Fig|31 The 
question we have to address now is whether it is possible 
to choose amplitudes A{^L,^fj,g) so that the function 
a{^L,^il, g\h) is equal to unity. 

To simplify notations, let us denote triples {^l,^r, g) 
by a single label s. The group G acts on these triples 
according to: 

s = ($L,$fl,.g) {h^Lh~\h<i>Rh-~\hgh^^) = h{s) 

(A4) 

Let us also introduce the function f{s\h) defined as: 



fis\h) 



( Xg1>i.(ft) Xj>fig-i(^) 



(A5) 
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The condition a(s\h) = 1 now reads: 
A{h{s)) 



A{s) 



fi-s\h) 



(A6) 



Now we show that solutions A{s) exist if and only if the 
function f{s\h) satisfies the following conditions: 

f{s\h) - 1 if h{s) = s (A7) 
fis\hih2)^f{h2is)\h,)f{s\h2) (A8) 

Note that the second condition is very similar to the fun- 
damental property (|17|) which enforces the existence of a 
local gauge symmetry group. In fact, using Ijl?!) . we can 
check that property (|A8() always holds. Property (|A7() 
is a new and quite strong condition imposed on gauge 
phase-factors X*(5); it is equivalent to the condition H45|) 
which in the main text has been derived as necessary. We 
now show that this condition is also sufficient. 
Let us consider the equation: 



Ajgis)) 
A{s) 



= fis\g) 



(A9) 



where f{s\g) and A(s) are complex numbers with unit 
modulus, s belongs to a set S on which the group G acts, 
and g Cz G. Here, the function f{s\g) is supposed to be 
known, and we are looking for "amplitudes" A{s). Let 
us denote by the stabilizor of s, namely it is com- 
posed of all the group elements h such that h{s) = s. 
From Eq. ljA9|l . we get the two conditions: 



(AlO) 
(All) 



f{s\h) = 1 if h{s) = s 
/(s|5i52) = f{g2{s)\gi)f{s\g2) 

The second equation is obtained from: 

Mgmjs)) ^ A{gig2{s)) A{g2{s)) 
Ais) Aig^is)) A{s) ' 

Note that these two conditions imply that: 
fislgh) = f{s\g) when h e Hs 



(A12) 



(A13) 



Let us now show that when these two conditions are 
satisfied, we can always reconstruct a system of ampli- 
tudes solving Eq. (jA9|) . From the form of this equation, 
we see that the various orbits in S for the action of G 
remain uncoupled. Let us then concentrate on one orbit 
generated by a fixed given element s. Elements of this 
orbit are in one to one correspondence with the left cosets 
gHs of G, since gh{s) = g{s) when h & Hg. Let us choose 
in each coset a representative element gn ■ Our problem is 
to find amplitudes A{gn{s)) knowing A(s). Let us define: 



A{gn{s)) =^(s)/(s|g„) 



(A14) 



From Eq. (|A13|) . we have A{g{s)) = A{s)f{s\g) for any 
g in G, since g may always be written as gnh with h G 



Hs- Now Eq. ljA9|l is an immediate consequence of the 
condition (|A11() . 

To conclude, when conditions ljA10|) and ljAll|) are sat- 
isfied, solutions of Eq. HA9|I can be constructed indepen- 
dently on each orbit of S for the action of G. On a given 
orbit, the solution is unique, up to the choice of A(s) for 
one particular element of this orbit. 



APPENDIX B: STRUCTURE OF THE GROUND 
STATE ON A TORUS 



The distinguishing feature of a torus is the appearance 
of non-trivial closed loops that are classified according to 
their winding numbers associated with two fundamen- 
tal cycles 72; and 7^, chosen with a common origin O. 
As a consequence, even if local fiuxes vanish on all pla- 
quettes, the fluxes and associated to and 
may not vanish. These global degrees of freedom are 
the origin of the topological degeneracies exhibited by a 
large class of gauge-invariant models on a closed space 
with non-trivial topology^^. How are these degenera- 
cies affected by the presence of a Chern-Simons term? 
Precise formulas for the Hilbert-space dimension of the 
pure Chern-Simons theory with a finite gauge group G 
have been given^" in terms of the action via functions 
a{h,k,l) in H^{G,U{1)). In this paper we consider a 
much larger Hilbert space induced by all allowed clas- 
sical gauge configurations, \{gij}). We may expect to 
recover the pure Chern-Simons Hilbert-space (which di- 
mension is independent of the system size) by projecting 
the full Hilbert-space onto the gauge-invariant ground- 
state of some local gauge-invariant Hamiltonian. This 
goes beyond the scope of the present paper because our 
goal here is classification and basic properties of such 
local gauge invariant Hamiltonian. However, it is in- 
structive to understand how, in the present formalism, 
a Chern-Simons term affects the gauge dynamics in the 
topological fluxless sector on a torus. 



FIG. 6: (Color online) For a square lattice on a torus, funda- 
mental oriented cycles 7^, and 7^ are depicted as black dashed 
lines. The vertical string of horizontal blue lines represents 
links carrying the group element $2, . Similarly, the horizontal 
string of vertical red lines represents links carrying ^y. On 
all other links ij on the square lattice gij = e, so we have not 
represented them not to overcrowd the figure. 
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♦ 





FIG. 7: (Color online) The upper line illustrates intermediate steps in the action of Cy{gy) on the template The 
plain (double) vertical red lines carry the element $y (gy^y) whereas the plain horizontal blue lines carry Successive 
locations of the fluxon moving from left to right are represented by a black cross. The lower line describes the action of 
the same string operator on the new state Cx{gx)\^x,^y) and the double horizontal blue lines carry now Qx^x- One goes 
from the upper to the lower line by a gauge transformation attached to the sites indicated by a green circle. Since this 
transformation commutes with the fluxon motion operators, these processes generate a contribution to X{^x, ^y \gx,gy) defined 
in section |n| equal to the product of the phase factors due to the two gauge transformations shown by thick black arrows. This 

factor is equal to \ x -i (ffa:)x<E.-i -i* {g^)) ■ ^ second contribution to M^x,^y\gx,gy) comes from the comparison between 
the actions of Cx{gx) on \$x,^y) and Cy{gy)\^x,^y}- This process is not shown on the figure, and yields a term equal to 

(xsx {gv)x^-\^^Jgy)) 



Starting from the trivial configuration gij = e, we may 
induce a new fluxless state in the bulk but with a non- 
trivial along by creating a string of parallel horizon- 
tal links carrying each the group element as shown on 
Fig. IHl This string may be viewed as the result of a suc- 
cession of elementary processes discussed in section IIVI 
first the creation of a fiuxon pair, then the mo- 

tion of one of them along and finally annihilation of 
the two fluxons. As we have seen, each of these processes 
may be described by a gauge-invariant operator. Let us 
denote by Cx{^x) the product of all these operators in- 
volved in the creation of the string. Similarly, we define 
Cy{^y). More generally, we may start from a state al- 
ready characterized by a pair of fluxes {^x,^y)- Since 
and 7y commute up to homotopy, we require that $^ 
and $y commute in G. The string picture allows us to 
construct operators Cx{gx) and Cy{gy) such that: 



($„$,) ''^^gA,*,) 
($.,$,) '^'($,,3,$,) 



Denoting by \^XT^y) any state with the pair of fluxes 
($a;, ^y), we show below that Cx(gx) and Cy{gy) do not 
commute, but their actions on a state are related 

by phase-factors: 



Cy{gv)Cx{gx)\^x.^y) = \{^x,^y\9x,9y) 

X Cx{gx)Cy{gy)\^x,'^y) 



\{^x,^y\9x.9y) 



Xgx {9yfx,^-^g^,i.S9yf 

Xg-^{9x?X^-^g-^^^{9x) 



Note that in order for all four states \^x,^y), \9x^xT^y) , 
\^x,9y^y)^ \9x'^x,gy'^y) to bc deflucd, the following con- 
straints have to be imposed: 



'^x^y 


= ^y'^x 


(B5) 


'^x9y 


= 9y'^x 


(B6) 


9x^y 


= '^y9x 


(B7) 


9x9y 


= 9y9x 


(B8) 



This relation ljB4p is reminiscent of the phase factors in- 
volved in permuting electrical field operators attached to 
adjacent links for Abelian models with a Chern-Simons 
terr 



19.27 



(Bl) Let us now briefly explain how to derive the rela- 

tions l|B3p and (|B4|) expressing the commutation rules be- 
(B2) tween Cx{gx) and Cy{gy). The first step is to choose tem- 
plate states for a pair of fluxes {^x^ ^y) along elementary 
cycles 'jx and 7y. These are depicted on Fig. El Let us 
now compare the action of Cx{gx)Cy{gy) and Cy{gy)Cx{gx) 
on this template \^x,^y)- The situation is illustrated on 
Fig. [71 Note that we do not need explicitely the values of 
the fluxon hopping amplitudes <i>j^, g), but simply 

the fact they are local and that fluxon moving opera- 
CBS) tors commute with gauge transformation. The value of 
X{^x,^y\9x,9y) obtained in this way (see Eq. (jB4p ) is 
(B4) valid not only for our template state \^x, but for any 
state deduced from it by a gauge transformation. 
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